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Abstract. We determine the spectral curve of charge 3 BPS su{2) monopoles with C3 
cyclic symmetry. The symmetry means that the genus 4 spectral curve covers a (Toda) 
spectral curve of genus 2. A well adapted homology basis is presented enabling the theta 
functions and monopole data of the genus 4 curve to be given in terms of genus 2 data. 
The Richelot correspondence, a generalization of the arithmetic mean, is used to solve 
for this genus 2 curve. Results of other approaches are compared. 
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1. Introduction 

The first order Bogomolny equations, 

1 3 

(1.1) ^» = 2 5Z '^i'^^'" = ^^^^ 

are rather ubiquitous. They arose while studying a limit of Yang-Mills-Higgs gauge theory 
in three space dimensions in which the the Higgs potential is removed but a remnant of 
this remains in the boundary conditions associated with ( 1.1 1. Here F^j is the field strength 



Date: June 25, 2010. 



2 



H. W. BRADEN, ANTONELLA D'AVANZO, AND V. Z. ENOLSKI 



associated to a gauge field A, and <i> is the Higgs field. The same equations may also be 
viewed as a dimensional reduction of the four dimensional self-dual equations upon setting 
all functions independent of and identifying ^ = A4] they are also encountered in su- 
persymmetric theories when requiring certain field configurations to preserve some fraction 
of supersymmetry. Just as the self-duality equations admit instanton solutions in four di- 
mensions, the Bogomolny equations possess topological soliton solutions with particle-like 
properties, known as magnetic monopoles, and these have been the subject of considerable 
interest over the years |MS04j . Early on a curve C was found to be associated to these BPS 
monopoles. Indeed the same curve arose by (at least) two different routes both with origins 
in instanton theory. Whilst considering the Atiyah-Ward instanton ansatz in the monopole 
setting Corrigan and Goddard |CG81) encountered C and this was given a twistorial de- 
scription by Hitchin fHit82]. Just as Ward's twistor transform relates instanton solutions on 
to certain holomorphic vector bundles over the twistor space CP'^, Hitchin showed that 
the dimensional reduction leading to BPS monopoles could be made at the twistor level as 
well and the curve C naturally lies in mini-twistor space, C C TP^. The second appearance 
of the curve C is closely connected with integrable systems. Nahm gave a transform of the 
ADHM instanton construction to produce BPS monopoles |Nah82j and the resulting Nahm's 
equations have Lax form with corresponding spectral curve C. Hitchin |Hit83) proved that 
all monopoles could be obtained by Nahm's approach provided the curve C was subject to 
certain nonsingularity conditions. Bringing methods from integrable systems to bear upon 
the construction of solutions to Nahm's equations for the gauge group SU{2) Ercolani and 
Sinha |ES89| later showed how one could solve (a gauge transform of) the Nahm equations 
in terms of a Baker- Akhiezer function for the curve C. Thus given a curve C the machinery 
of integrable systems allows one (in principle) to construct solutions to Nahm's equations 
and thence monopoles |BE09Aj . 

The problem in the approach just described, and to which this paper is devoted, is in 
constructing the curve C: some of the conditions necessary for the regularity of solutions just 
alluded to impose transcendental constraints on C, and we presently lack analytic means for 
solving these. One such constraint comes about by requiring the periods of a meromorphic 
differential on C to be specified. This type of constraint arises in many other settings as well, 
for example when specifying the filling fractions of a curve in the AdS/CFT correspondence 
[KMMZj . finding closed geodesies on an ellipsoid |AF06j or constructing harmonic maps 
— >■ [Hit90| . The second type of constraint is that the linear Sow on the Jacobian of 
C corresponding to the integrable motion only intersects the theta divisor in a prescribed 
manner; equivalently this may be expressed as the vanishing of a real one parameter family of 
cohomologies of certain line bundles on C. While techniques exist that count the number of 
intersections of a complex line with the theta divisor we are unaware of anything comparable 
in the real setting :BE09Bj . Thus the application of integrable systems techniques to the 
construction of monopoles (and indeed more generally) encounters two types of problem 
that each merit further study. 

In the present paper we will simplify then solve these constraints by imposing spatial 
symmetries on the monopole. Imposing symmetry reduces the number of constraints to be 
solved for and here we shall focus on charge 3 monopoles with cyclic symmetry C3. The 
spectral curve of a charge n monopole may be expressed in the form 

Piv, C) rf + r7""iai(C) + • ■ • + V^'dnMO + ■■■ + '7«n-i(C) + ^^(C) = 0, 

where ar(C) (for 1 < r < n) is a polynomial in of maximum degree 2r. Hitchin's construc- 
tion involves three constraints on the curve. The first (HI) requires the curve C to be real 
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with respect to the standard real structure on T¥^, 
(1-2) t:{C,v)^{-^,-^), 

the anti-holomorphic involution defined by reversing the orientation of the lines in M'^. As 

a consequence the coefficients of the curve satisfy ar{C) = (— l)'^C^'^ar(— 1/C)- Imposing 
spatial symmetries via fractional linear transformations of TP^ simplifies the curve. Long 
ago monopoles of charge n with cyclic symmetry C„ were shown to exist |OR82j and these 
correspond to curves invariant under (?7, C) ~^ i^li'^Q) (where lj — exp(2i7r/n)). Imposing 
Hitchin's reality conditions and centering the monopole (setting ai — 0) then gives us the 
C„ symmetric spectral curve in the form 

r;" + a2??""'C' + • . ■ + a«C" + /?C'" + (-1)"^ - 0, a, G R. 

By an overall rotation we may choose /3 real and so the charge 3 spectral curves C we will 
focus on in this paper have the form 

(1.3) + ar/C^ + /^C'' + 7C^ - /3 = 0. 

The remaining two constraints of Hitchin on this curve are the transcendental constraints re- 
ferred to above. The first of these (H2) may be expressed in the following manner jHMROOj : 
given a canonical homology basis {a^, b^} for the curve C there exists a 1-cycle es = n-a+m-h 
such that for every holomorphic differential 

(1.4) „. ^.<-' + ft(cg-^+,.. + ft-,(o ,^_ 

Cs 

Dually the vector 

U = -n H — f m 
2 2 

(where f is the period matrix of C) is a half-period. The vector U is known as the Ercolani- 
Sinha vector [ES89 and may be expressed as the periods of a meromorphic differential. 
These conditions, known as the Ercolani-Sinha constraints, impose (n — 1)^ (the genus of 
C) transcendental constraints on the curve. The remaining constraint (H3) is that for a 
special vector K in the Jacobian the linear flow XU — K intersects the theta divisor only 
at A = and 2. The consequences of assuming symmetry is that the spectral curve covers 
another curve, the quotient curve by the symmetry. In the case of cyclic symmetry we have 
an n-fold unbranched cover tt : C — >■ C of a hyperelliptic curve of genus n — 1 which is the 
spectral curve of the a„ affine Toda system, and |BralO| shows how both the constraints 
(H2,3) reduce to become constraints on the reduced curve. In particular 

XU-K = TT*{XU - Koo+ + e). 

where U expresses the Ercolani-Sinha constraints of the curve C and the other quantities 
will be defined later. Thus in this paper we have a 3-fold unbranched cover of the of the 
curve C given by 

(1.5) ^ {x^ + ax + "f)^ 

To solve the remaining transcendental constraints for the reduced curve C we now use 
several pieces of research. First the work of [BE061 |BE09BQ identifies the solutions of 
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for the class of curves (1.3) with a = and |BE09Bj shows that the only solutions of the 



Hitchin constraints are for the curves 

1 nmi) 



(1.6) + ;t(C6 ± 5^2(3 _ 1) ^ 0, 



6 25 7r5 



These correspond to tetrahedrally symmetric monopoles and for each sign of ±5\/2 there is 
a unique vector U. Next we use the work of [ HMM95] . Here cyclically symmetric monopoles 
(and more generally, those with Platonic spatial symmetries) were reconsidered from a va- 
riety of perspectives. Cyclically symmetric monopoles form 4-dimensional totally geodesic 
submanifolds A^J^ of the full moduli space of charge n monopoles, where < I < n. Ignor- 
ing the rotational degrees of freedom these then yield one dimensional submanifolds. By 
considering the rational map description of these monopoles Hitchin, Manton and Murray 
were able to further specify the Ai^ which may be viewed as orbits of geodesic monopole 



scattering. For charge 3 there were five loci of spectral curves of the form (1.3 1. Of these 
loci, four were isomorphic: at one end asymptotically one has — 270^ (with /3 of either 
sign) and 7 = while at the other end a = 7r^/4 — 36^, /? = and 7 = 26(6^ -f 7r^/4) (with 
h of either sign). Half-way along this is the tetrahedrally symmetric monopole, the four loci 
corresponding to four distinct orientations of the tetrahedron. The final locus corresponds 
to the family of curves with 7 = and where the symmetry is enlarged to the dihedral 
symmetry D3: asymptotically we have = 27/?^ (with /3 large and positive at one end and 
negative at the other) and half-way along this there is the axisymmetric monopole. We use 
this work as follows. Because the Ercolani-Sinha vector U is discrete, this will be constant 
for each of the loci emanating from the tetrahedrally symmetric points. Starting then at a 
point (a,/3,7) corresponding to a tetrahedrally symmetric monopole we deform away from 
this by solving the (reduced) Ercolani-Sinha constraint for the given fixed U . Thus we will 
obtain the loci 7 7^ 0. This idea is similar to that used by Sutcliffe |Sut97j when obtain- 



ing numerical approximations to (1.3) by analysing the Nahm equations and using the fact 
that the tetrahedral monopole was on one of the loci. In deforming from the tetrahedrally 
symmetric points we will use a genus 2 variant of the arithmetic-geometric mean (AGM) 
(that will be described more fully in the sequel). Although this deformation is defined by 
the (reduced) constraint (H2) the ensuing loci must also satisfy (H3) for dimensional rea- 



sons. Thus we arrive at the spectral curves (1.3) with 7 7^ that describe C3 symmetric 



monopoles. We remark that the 7 = spectral curves further cover an elliptic curve and 
these are amenable to a different analysis that will be given elsewhere. 



An outline of the paper is as follows. In section 2 we study the curves (1.3 1.5) in some 
detail determining those quantities needed to reconstruct the Baker- Akhiezer functions of 
the integrable systems approach. Critical here is determining an homology basis that reflects 
well the symmetries of the curve C. Such a basis both relates and simplifies the forms of the 
period matrices of both C and C; it also reduces the numbers of periods to be calculated to 
construct the full period matrices. Perhaps the nicest feature of this homology basis (and 
that induced on C) is that it enables us to make use of a remarkable factorisation theorem 
due to Accola and Fay |Acc71[p^ay73| and also observed by Mumford. This allows the theta 
functions of C to be described in terms of the theta functions of C for the parts of the Jacobian 
that are relevant for us |BralO| . By the end of section 2 we have reduced the construction 
of cyclically invariant monopoles (with 7 7^ 0) to questions about a genus two hyperelliptic 
curve C. Section 3 then discusses the restrictions the Ercolani-Sinha constraints place on C. 
The Ercolani-Sinha constraints are shown to reduce to the single constraint j^dX/Y — 
on a scaled form of C, = [X'^ + + 5)^ + 4 (where (0,5) := (a//3^/^, 7//?)), and here 
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c :— 7r(es) is known. Thus the problem has become one of understanding the periods of this 
(scaled) genus two curve as a function of (a,g) with the Ercolani-Sinha yielding g = g{a). 
We will solve this transcendental constraint numerically using a genus two variant of the 
arithmetic-geometric mean due to Richelot. Section 4 recalls this theory and describes an 
extension needed for the curves relevant here, which have complex conjugate branchpoints. 
Section 5 then implements this. We conclude with a discussion. 

2. The curve 



In this section we consider the curve C ( 1.3 1 and the quotient curve C in more detail. After 
describing the curves we shall construct homology bases that enables several simplifications. 
In particular both the period matrices and vectors of Riemann constants will be described 
for these bases. Throughout we will set p — exp(2j7r/3). 



2.1. Branchpoints and monodromy. The curve ( |1.3[ ) has genus 4 and is not hyperelliptic. 
A basis for the holomorphic differentials may be taken to be 



(2.1) Ui — 2 I A2 ' "2 — „ 2 I 72^ "3 — 72 ' — „ 2 i ^2 ' 

Our curve C may be viewed as a 3 sheeted cover of with 12 ramification points, whose 
('-coordinates are given by 



B 



1,2 



- — f-97 - 2z^/3a3/2 ± A^' 

8/3 V ' ^ 



18/3 



(2.2) Bl, = ^ (-97 + 2^V3ay' ± A'J' 

%4fe-/(S^)'/', J-l,...,4, fc = 0,l,2, 



where 

A± = 324/3^ - 3(2a^/2 ± SiVSjf. 

In Figure [l] we give a qualitative sketch of the branchpoints; the general properties of C 
do not change with the parameters unless a = 0, which is a degenerate case that will be 
examined separately. The monodromy around each branch point is found to be 

Bi, i?2, By, Bg, — > [1,3], i?3, i?4, Bg, Biq — > [1,2], i?5, Bq, Bh, B12 — > [2,3]. 

In addition to the real involution ( |1.2[ ) and cyclic symmetry 
(2.3) a: ^ {pC PV) 

the curve possesses the inversion symmetry 



(2.4) </.:(C,r?) 



1 _^ 



The branchpoints (2.2) form four orbits under the cyclic symmetry according to 

i>l > i'5 > ^9, -02 > -D6 > -DlO, 

-D3 > ^7 > i"!!, -fc>4 > i>8 > i'12- 

The case a = 0. The case a = is the curve studied in |BE06| . The corresponding Riemann 
surface also has genus 4, but now only six branchpoints (i = 1 . . . 6), 
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[1 ,2] 



,^/\/\/\/\/ [1 ,3] 



[2,3] 



Figure 1. Branchpoints and monodromy for C 



A2 = pAi, A3 = p^Ai, 



A5 — pA4, Ag — p^Xa 



where A ~ 27(12/3^ + 7^). Indeed, letting a — > in (2.2), we sec that the branchpoints Bi 
coUide pairwise to give the \i (see Figure [2]). In Figure [3| we also give the monodromy which 
is [1, 2, 3], the same for every branchpoint; this can be seen from an explicit calculation, but 
also by taking the limit of the monodromies of C. 



2.2. The quotient with respect to C3. We may form the quotient curve C — C/a with 
the covering map 

1 



(2.5) 



C 



and the curve C is given by (1.5). It is a genus 2 (hence hyperelliptic) Riemann surface and 
TT is an unbranched covering. Viewing C as a 2-sheeted cover of the Riemann sphere, it has 
six branchpoints whose x-coordinates are 

2p2a 



(2.6) 



where 



1 

- no- 



2a 



Bi 



1 



SI 



2p^a 

2 ' 

2pa 



Bi 



Bi 



6 - 



2a 
SI 



1 



2 A3 



■p^S 



2pa 
61 



S± =-1087-216/3^ + 12^/12^3 + 81 (7±2/3i)^ 
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3 -] 



2 - 




-2 - 



^ Branchpoints >^ ^ Branchpoints B ^ varying a 



Figure 2. Branchpoints for a -> 




Figure 3. Branchpoints and monodromy for a — > 0. 

As (5_ = (5+, these branchpoints can be spht into complex conjugate pairs 
(2.7) B^^Bf, B^^Bl, Bl^Bl. 

The branchpoints Bi are not the images of the branchpoints of C under tt. Figure [4] shows 
again a quahtative sketch of the branchpoints for the curve and their monodromy, with the 
same choice of parameters of Figure [l] Because the branch points are not images of those of 
C we observe that there is httle difference in the quotient curves for a equalhng or differing 
from zero. 
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Figure 4. Branchpoints and monodromy for C 



A standard basis for the holomorphic differentials on C is given by 
(2.8) 



da; X dx 

Ui = , U2 . 

y y 



The differentials U2 and U4 on C are invariant under a and hence descend to differentials on 
C. In fact, 



(2.9) 



TT* Ui = -3U2, 



TT* U2 = — 3U4. 



This observation allows us to considerably simplify some integrals, and hence the period 
matrix. 

2.3. Homology bases. By choosing an homology basis well adapted to the symmetry at 
hand we may simplify many things. The aim of this subsection is to construct one such 
basis. For the case of our unbranched cover tt : C C it is known |Fay73| that there exists 
a basis {og, bo, Oi, bi, 02, b2, a^, b^} of homology cycles for C and {oo, bo, Oi, bi, } for C such 
that (for 1 < j < 3, < s < 3) 



(2.10) a^(ao) ^ ao, 

(2.11) 7r(ao) = ao, 



(J^cij) = hj+s, cr^(bo) = bo, c^''(bj) = b^+s 

n{aj+s) a-,-, 7r(bo) = 3 bo, 7r(bj+s) = bj. 



Here cr*(oo) is homologous to ao and indices are understood to be modulo 3. These require- 
ments do not uniquely determine a homology basis and we may impose some extra conditions 
on this basis. The condition we choose (and will see is possible) is that the limit a — !• of 
some cycles are mapped to some of the homology basis in |BE06j . This will enable us to 
relate the present work with |BE06j whose results we generalise. We also observe that given 
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two cycles di, bi with canonical pairing then we may simply define £2,3 via di+s — cr'*(di) 
(and similarly for 62,3)- Thus we seek two cycles di, bi such that 

d? lim di = af b? := lim bi = hf^''. 

a— ^0 a— 

Such cycles and their corresponding images under a are shown in Figure [5] alongside the 
a = limit of these. One further finds that we may choose do such that 

d° lim do = afEoe 



and (2.101 is satisfied. 



It remains to find the cycle bo and this is the most difficult. From (2.10) we see we wish 
a cycle invariant under a and having canonical intersections with the other cycles. Such is 
shown in Figure [6] alongside do and their a = limits. We record these results as 



Theorem 1. The homology cycles given in figures^ and^ are canonical, satisfy (2.10), 
and have smooth limit a — > 0. 

These cycles can be expanded in terms of "basic arcs" as follows. Denote the arc between 
the branchpoints Bi and Bj on sheet k by 

lk{i,i) a.TCk{B,,Bj), z 7^ j = 1, . . . , 12. 
Then we have the following 

di = 71(1, 2) + 72(2, 1), bi = 7i(3, 1) + 72(1, 12) + 73(12, 3), 
d2 = 72(5, 6) + 73(6, 5), b2 = 72(7, 5) + 73(5, 4) + 71 (4, 7), 

(2.12) da -73(9, 10) +71(10, 9), 63 = 73(11, 9) +71 (9, 8) +72(8, 11), 
do = 7i(3, 10) + 73(10, 9) + 7i(9, 8) + 72(8, 12) + 73(12, 3), 

bo = 7i(2, 8) + 72(8, 11) + 73(11, 4) + 71 (4, 7) + 72(7, 12) + 73(12, 2). 

and, as a — )• 0, 

d? = 71(1, 2) + 72(2, 1), b° = 71 (2, 1) + 73(1, 2), 
a" = 72(3, 4) + 73(4, 3), b° = 72(4, 3) + 71 (3, 4), 

(2.13) d° = 73(5,6) +71(6,5), b°= 73(6, 5) +72(5, 6), 
d° = 73(1, 2) + 71(2, 6) + 73(6, 5) + 72(5, 1), 

bSJ = 73(1, 2) + 71(2, 5) + 72(5, 6) + 73(6, 3) + 71(3, 4) + 72(4, 1). 



We may complete the specification ( 2.11[ ) of homology bases by projecting the cycles of 



Figures [5| [6| The fact that the branchpoints of C do not get mapped by tt to branchpoints 
of C makes the projection less straightforwarcj^ The results are shown in Figure [t] We 
therefore have a homology basis for the hyperelliptic curve C differing from standard ones. 
As we shall see however, the bases chosen allow us to simply relate the period matrices and 
other quantities of C and C. With the same notation as above, the arc expansion for these 

"^This has been implemented in Maple. 
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Figure 5. Cyclic homology basis 




Figure 6. Cyclic homology basis 
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Figure 7. Cyclic homology basis for the quotient curve C 



cycles is then 

ai =7i(2, 6) +72(6,2), 
bi =71(6,4) +72(4,6), 

( 2 14 1 

^' ' ao = 71(3, 4) +72(4, 6) +71(6,1) +72(1, 6) +71(6, 4) +72(4, 3), 

bo = 71 (3, 4) + 72 (4, 5) + 71 (5, 4) + 72 (4, 3) . 
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2.4. Period matrices. We shall now relate the period matrices of C and C and then use 
the symmetries of C to further restrict the periods involved. If {a^, b^} are the canonical 
homology basis introduced earlier and {uj} any basis of holomorphic differentials for our 
Riemann surface C we have the matrix of periods 



(2.15) 




A 



with f — BA^^ the period matrix. The period matrix and periods A, B are our focus here. 

To understand the connection between the period matrices f of C and r of C we may first 
focus on the a-normalized differentials. If are the d-nornialized differentials for C, then 



and wc find that 

(2.16) (a^)*vo=vo, (a^)*v, =v,; 

If V; are the normalized differentials for C, then 



'r*(v.) 



7r*(vi) =v, + (ct)*v, + ((72)*V, 



shows that 
and similarly that 

7I'*(Vo)=Vo- 

From (2.16) we have an action of C3 on Jac(C) which lifts to an automorphism of by 
(2.17) a{z) = a-(zo, zi, Z2, ^3) = (^0, ^3, 22, zi) 

With the choices above (things are different for b-normalization) we may lift the map tt* : 
Jac(C) ^ Jac(C) to tt* : ^ C^ 

7r*(z) = 7r*(zo,2:i) = {3zo,Zi,zi,zi) = z. 

If we denote points of the Jacobian in characteristic notation by 



(a, /3 e ti^gj^ 





ax 




ao 


ai 


ai ai 




Pk 


T 


3/3o 


/3i 


Pi Pi 



(2.18) 



2d 



The period matrices for the two curves are related by |Fay73| 

( a h b b\ 
b c d d 
b d c d 
\ b d d c J 

The point to note is that although the period matrix for C involves integrations of differentials 
that do not reduce to hyperelliptic integrals, the combination of terms appearing in the 
reduction can be expressed in terms of hyperelliptic integrals. This is a definite simplification. 
Further the Q function defined by f has the symmetries 

e(i|f) = e(cT^(z)|f) 
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for all i e C^. In particular, the Q divisor is fixed under C3. 

Now we turn to the symmetries of C to simplify the calculation of periods. If ip is any 
automorphism of C then ■0 acts on i/i(C,Z) and the holoniorphic differentials by 



V'* 



where 



e Sp{&,Z) and L e GL(4,C). Then from 



we obtain 
(2.19) 

Thus, for example, from the definition of our homology basis, 



A B\ (A 
C D [b 



L. 



( 1 




















\ 








1 


























1 

















1 
































1 





























1 


























1 


V 














1 





/ 



while a acts on the differentials (|2.1|) as follows 
(2.20) 



(7 Ul = 



CT U2 = U2, 



cr U3 p U3, 



cr U4 = U4. 



Let us denote the hi integrals of Ui by Zi and the corresponding hi integrals by Zi, and simi- 
larly those of U2, U3 and U4 by Xi, Xi, Wi, Wi and yi, Yi respectively. Then A = (z, x, w, y) 
(and analogously for B). A symmetry relates these various periods. The symmetry a re- 
stricts the matrices of periods to take the following form 

/ 



(2.21) 



A 





zi 



Xi 
Xi 
Xi 





Wi 

pwi 



For instance, we get 



zo = 0) Ul 



yi 
yi 



Ul 



B 





^2 ' 



Xi 







Yo 
Yi 
Yi 
Yi 



a Ul : 



p^Ui = P^Zq, 



and so zq = 0, while 



Z2 = f Ul 

'02 



Ul = 3) (7 Ul 



p^Ui = Zi, 



SO leading to (2.211 



We now use the other symmetries to further restrict the matrices of periods. The chief 
difficulty in this approach is in calculating the actions on the homology. In the present 
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setting we find that 
(2.22) 

/ 



:= 



with 

(2.23) 

(2.24) 



1 




















\ 








-1 




















-1 





























-1 


























-1 





























-1 

















- 


-1 





























-1 / 




Ui = 










-U2, 




T* 


Ui = 


U3, 




T U2 




-U2, 





/ 2 
10 11 
110 1 
1110 
6 111 
10 
10 
\ 1 





-1 
-1 
-1 
-2 






-1 

1 
1 

-1 

-1 
-1 



-1 \ 
1 
1 


-1 
-1 
-1 

/ 



if Ug 
T U3 



Ul, 
Ul, 



if U4 
T*U4 



-U4, 
-U4. 



The if symmetry simphfies the matrices of periods to the form 



/ 



A 



xo 2/0 \ 

xi -p^zi yi 

xi -Zl yi 

\ pzi xi -pzi yi ) 

while the real involution relates the entries via 



B 







Xi 
Xi 
Xi 





-Zl 
-pZi 



Yi 
Yi 



pZl - Zl, 



Wi = -(zi 
3Yi = 2yo + yo 
Yq ^ 2yi + yo 



Wl) 



P Zl - Zl, 



Vi+IYi- 



(2.25) Zi = -{wi+zi) 

(2.26) 3X1=2x0 + ^0, 

(2.27) Xq = 2x1 +xa + xi+ 2Xi , 

Thus all the periods are determined in terms of zi, xi, xq, yi and y^. Finally there is the 
bilinear relation 

= xqYo - yoXo + 3{xiYi - yiXi). 
Calculating the period matrix f = BA^^ we obtain the form (2.18) with 



(2.28) 



2Zi 

3 Zl 



Xpyi 
IxoYi 



xiYi - Xiyi 



Xiyo 



xiya - 
_lZi 
3 Zl 



IxoYi 



Xiyo 



3 Xiyo - xoVi 3 Zl 3 a;iyo - x^yi 

2.4.1. Weierstrass-Poincare reduction. We remark in passing that the symplectic matrix 



(2.29) 



T ^ 






1 


1 


1 














1 






































-2 


1 


1 








1 


-1 





























1 




















1 




















-1 





























-1 


1 






A 
C 



B 
D 



/ 



transforms the period matrix (2.181 as 



f ^ {C + Dt){A + Bt)- 



Q 
1 
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where Q = Diag(— 1/3, 0) and 



c—d 
6 



2(c-d) 



From this we deduce that 

(2.30) Im(c - d) 7^ 0. 

2.4.2. The antiholomorphic involution for C. We have seen that the spectral curve C has 
real structure ( 1.2 1. This real structure is inherited by C where we have the antiholomorphic 
invohition 



(2.31) 



/ : {x,y) {x,-y). 



The effect of this is on the homology above is to reflect in the x-axis and change sheet. 
Specifically we find that 



(2.32) 



Mr' 



/2 -3\ 

12-32 
2 1-2-1 

yi -2 J 



M^, = Id, Mr'JM^, 



and where the last identity reflects that Mt' is antiholomorphic. On the holomorphic dif- 
ferentials we have the simple action 



T Ul 



-ui. 



T U2 



-U2. 



2.5. The Fay-Accola theorem. Having established homology bases (2.10 2.111 and the 
relationship these entail for the corresponding period matrices (2.18) of the curves, we next 
recall the striking theorem of Fay and Accola applied to our present setting. 

Theorem 2 (Fay-Accola). With respect to the ordered canonical homology bases {di,bi} 
constructed above and for arbitrary z =G 

(2.33) 



0[e] 






"0 0' 

- 


(^;t) 



we have that 

co(t) 



is a non-zero modular constant co(t) independent of z. Here f and r are the a-normalized 
period matrices for the respective curves given in the above bases and 



7r*(e) 



1 



2-3 



,0 -(1,0,0,0) 



0. 



The significance of this theorem is that for flows on the Jacobian of C that arise as 
pullbacks of flows on the Jacobian of C we may reduce the theta functions to those of the 
hyperelliptic spectral curve. We have stated in the introduction that such a connection 
holds, 

XU -K^ Ti*{XU - Koo+ + e), 
and we now describe the quantities appearing in this. 



ON CHARGE-3 CYCLIC MONOPOLES 



17 



2.6. The vector of Riemann constants. To construct the Baker-Akhiezer function for 
monopoles there is a distinguished point K e Jac(C) that Hitchin uses to identify degree 
g — 1 hne bundles with Jac(C). For n > 3 this point is a singular point of the theta divisor, 
K e Osinguiar [BEOBj . If wc dcuote the Abel map by 



MP) 

then 



p 

Ui 

Q 



(2.34) K ^K^+An\{n-2)y^6Dk 



k=l 



Here Kq is the vector of Riemann constants for the curve C and obk are the points above 
infinity for the curve. If IC^ is the canonical divisor of the curve then Aq{ICq) = —2Kq. 



The righthand side of (2.34) is in fact independent of the base point Q in its definition. Let 



7r(6bfe) = 0O+ € C denote the projection of the points at infinity. Then |BralO| shows that 
(2.35) K ^Tr*{K^^)~e = Tr*iK^^-e), 

where the half-period e has been identified in ( |Fay73| ) . Thus we need to calculate the vector 



of Riemann constants (for the homology bases constructed) for the genus 2 curve (1.5) and 
where the basepoint for the Abel map is oo+. It will be easier for our calculations to choose 
one of the branchpoints, say Bi, and then to obtain the vector of Riemann constants with 
respect to oo+ by using the relation 

(2.36) K^.^=AaoABi)+KB„ 
where Aqo^ is the Abel map with basepoint oo^. 

2.6.1. The vector Kbi- We begin by expressing the integrals over our homology cycles in 
a simple form as integrals between branch points. Let ji{j,k) denote the path going from 
branchpoint Bj to Bk on the cut plane of Figure |4] corresponding to sheet i. With this 
notation the cycles ( 2.14[ ) can be expressed as 



ai = 71 (2, 1) + 7i(l, 6) + 72(6, 1) + 72(1, 2), 
bi = 71(6, 5) + 71(5, 4) + 72(4, 5) + 72(5, 6), 

ao = 71 (3, 4) + 72(4, 5) + 72(5, 6) + 71(6, 1) + 72(1, 6) + 71(6, 5) + 71 (5,4) + 72(4, 3), 
bo = 71 (3, 4) + 72(4, 5) + 71(5, 4) + 72(4, 3). 

These expressions may be further simplified using the hyperelliptic involution J : (x, y) — > 
(x, —J/), giving 7A;+i(2j, 2j + 1) = Jjk{2j, 2j + 1). Also, as the sum of the cycle cncirchng 
Bi and B^ together with the cycle encircling B^^ and ^4 and the cycle encircling B^ and Bq 
is homologically trivial, we find that relations 



UJ — UJ + UJ, 

7fc(6,5) Jjkil,^) J-YkiaA) 



for any holomorphic differential w. Similar expressions result from other homologically trivial 
choices of cycles. These yield 
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UJ = 2 



W = 2 



7i(l,2) "'71(3,4) 



71(3,4) 



71(1, 2) 



7i(5,4) 



71(3,4) 



7i(5,4) ^71(3, 2) 



71(5,4) 



bo 



71(3,4) 



71(5,4) 



and these may be inverted to obtain the integrals between branchpoints. If the ui are taken 
to be a-normaUzed differentials we have 



1 



71(1,2) 



71 (4,5) 



71(2,3) 



71 (3,4) 



^ l(e(o)+e(i)+. 



= i(e(°)+eW), 



i(e(o)+e(i)+T(i)) 



71(5,6) 



where t^^^ and e*^*^ are the appropriate rows of the period and identity matrices. Thus one 
can easily deduce the image of each branchpoint under the Abel map (with basepoint Bi). 
We obtain (together with their characteristic form) 



^Bi(Bi) -0 



1 









1 

1 

1 



^B.(S2)-^(r(°)+r«) 



^B.(i?4) = ^(e«+r("))^ \ 
, ^5,(i?6)^i(eW+r(")+rW)^ \ 



1 1 



1 

1 

1 1 

1 



Following an argument of Farkas and Kra ([FK8O] VII. 1.2), the vector of Riemann constants 
takes the form 



(2.37) Kb, = -{AbAB^) + AbAB^)) = i(e(°) + + t«) ^ ^ 



1 

1 1 



2.6.2. The vector K^oj^ ■ Once we have calculated the vector of Riemann constants with Bi 
as basepoint, we can change its basepoint making use of equation (2.36). One finds 



(2.38) 

and consequently 
(2.39) K 



-4oo,(Si) = ^eW + V)^ 



1 

2 n 



00. - ^e(i) + -e(2) + 1^(1) + 1^(2) 
°°+ 6 2 2 2 



1 1 

2 2 
1 1 

6 2 



2.6.3. The case a = 0. In the case where a = 0, we have that 3 J G A for any branchpoint, 

Bi 

and moreover jBEOe] that Ag^ (j2L 

obi. I =0. Thus in this case we have 



\fc=i / 



^Bi = 7r*(Koo+) - e = TT*{Koo+) 
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which yields 



1 1 
1 1 



This coincides with the result of |BE06) derived by other methods. 

3. The Ercolani-Sinha conditions 

Here we shall express the transcendental Ercolani-Sinha constraints on the curve C as 
conditions on the curve C and then describe our strategy to solve them. 

With the ordering of the differentials (2.1 1 the Ercolani-Sinha conditions (1.4) take the 
form 

=-2(0,0,0,1), 



(n, m) 



where es = n • d + m • b. Now substituting (n, m) = (no, ni, n2, n-s, toq, 'tii, m2, 7713) directly 
into (2.21) and making use of (2.28 2.30) we may deduce that the Ercolani-Sinha vector 
takes the form 

(n, m) = (no, n, n, n, mo, m, m, m), 

and thus es is fixed under the spatial symmetry: a{ts) = es. (This result was obtained 
more generally via a different argument in [BralO| .) With this simplification we find the 
remaining equations encoded in the Ercolani-Sinha conditions take the form 



Xi 



yo 
yi 



{mo, 3m) 



-2(0,1) 



Now using (2.28) we have that 



-'^0 Yq 
Xi Yi 



Xl 



yo 
yi 



a/3 
b 



b 

c + 2d 



3 
1 



Upon noting (2.9) and that the periods of A and B were constructed from u* we obtain 
Theorem 3. The Ercolani-Sinha constraint on the curve C yields the constraint 



(3.1) 



(no, 3n, 3mo, 3m) 



B 



= -2(0,1) 



on the curve C with respect to the differentials —dx/{3y), —xdx/(3y) and the homology basis 
{oo,bo,ai,bi,}. 



^(y,n) + ^(mo,m)r. 



We remark also that we have 

U = TT*{U), U 

If we define the cycle 

(3.2) c := 7r(es) — noOo + 3nai + 3mobo + 3mbi 
then the Ercolani-Sinha constraints may be alternately expressed as 

(3.3) 6/3o 



j) n* (/3oU2 + /3iUi) j) (/3oU2 + /3iUi 



At this stage then we have reduced the Ercolani-Sinha constraints on the curve C to 
analogous conditions on the curve C. We now use the approach outlined in the introduction. 
We know from [BE06j the values of {n,m) of the Ercolani-Sinha vector for the curve (1.6) 
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for both signs. After changing from the homology basis of that work to that of the present 



paper we obtain for the two cases of (1.6), 



(3.4) 



(no, n, mo, to) 



(4,1,-3,1) 
(5,1,-3,0) 



+5^2, 
-5V2, 



and so we know the cycle c for each of the two loci associated to the tetrahedrally sym- 
metric monopoles. We also remark that just as T^,{^s) = — e s IHMROOj we also have that 
T^(c) = — c. Thus for the +5-\/2 values above and using (2.32) appropriate to this we 
have (4, 3, -9, 3)Af^. = -(4,3,-9,3) and similarly for -5V2 we have (5, 3, -9, 0)1/^- = 
-(5,3,-9,0). 

Now if we now make a change of variable 



y = l3Y, a 



then 



y2 = (x3 + aX + 5)2+4, 



dx 

Ul = — =13 
V 



■2/3 



a 
dX 

Y y 



7 
/3 



xdx _ XdX 



and the Ercolani-Sinha constraints take the form 

dX 

1^' 



(3.5) 
(3.6) 







6/3 



1/3 



XdX 
Y ' 



(We denote by c the cycle for both the scaled and unsealed curves.) The first of these 
equations may be viewed as defining g = g{a), and then for this solution the second gives 
us /? = /?(<?) • Thus solving the Ercolani-Sinha constraints has reduced to determining the 



relation between g and a given by ( 3.5 ). Thus to solve the Ercolani-Sinha constraints we need 
to be able to compute periods of these hyperelliptic integrals. We shall do this numerically 
using a variant of the arithmetic-geometric mean used to rapidly and accurately compute 
periods of elliptic integrals. We shall turn to this in the next section. 



4. The AGM method 

In this section we shall recall the connection of the arithmetic-geometric mean (AGM) 
to evaluating elliptic integrals and Richelot's generalisation of this to the genus two setting. 
This latter work has been most studied in the setting where the hyperellptic curve has real 
branch points and we shall need to extend this discussion to the case with pairs of complex 
conjugate branch points relevant to the monopole setting. 

4.1. AGM: the elliptic case. While the origin of the AGM method dates back to Lagrange 
it was Gauss who truly initiated its investigation. A large part of what is known today 
seems to be due (or at least known) to him (for historical notes see e.g. \Cox84\ ). Let a> b 
be positive real numbers. The arithmetic- geometric mean of these numbers, denoted 
M(a, b), is the common limit of the sequences defined as follows: 

oo = a, bo = b, 

(4.1) a„+i = " ^ ^" , bn+i = V anbn- 
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These two sequences satisfy 

ao > ai > . . . > an > a„+i > . . . > 5„+i >&„>...> 61 > 60, 
which ensures the existence of a common hmit 

Um a„ = hm 6„ = M(a,b). 

Indeed ^ 
whence 

< a„ - 6„ < 2-"(a- 6), 

which ensures rapid convergence (which is relevant in the present work). 

The remarkable observation of Gauss was the connection of the Arithmetic-Geometric 
Mean with elliptic integrals. 

Theorem 4. AGM |Gau99) Let a, 6 G M-i- and let M (a, b) he their arithmetic geometric 
mean, then 

d0 _ 
Jo cos2 + 62 sin^ ~ 2M {a, b) 

The theorem may be understood in terms of maps Gn ■ £n ~^ ^n+i between the elliptic 
curves 

(4.2) £"„ : yl = Xn{xn-al){xn-bl). 

all elliptic integrals of the form J^^ '^^^ expressed in terms of an 



Using Theorem 

appropriate aritnmetic geometric mean after various change of variables. We also remark 
that the restriction a, 6 e may be extended to a,b £ C* , a ^ ±6, with further discussion 
of the square roots taken in the geometric mean. We shall not need this extension here. 

4.2. Richelot and Humbert: the genus 2 case. Richelot |Ric36[ lRic36j extended to the 
hyperelliptic case Gauss' connection of the AGM with elliptic integrals. Humbert [HumOlj 
later gave another view of this re-interpreting Richelot 's findings in terms of the duplication 
formulae of 2-variable theta functions, i.e. isogenics (of type (2,2) ) on Abelian surfaces. 
We will follow here the modern exposition of Richelot 's work given by Bost and Mestre in 
|BM88| which describes Richelot 's "changes of coordinates" in terms of a correspondence 
(see below). 

At the outset we note that the Richelot-Humbert construction is only given for the case 
where the genus 2 curve (represented as a two-sheeted cover of with six branchpoints) 
has all real branchpoints. This manifests itself in what follows by using the ordering of the 
reals. Just as the AGM of two complex numbers is correspondingly more complicated than 
the real setting the implementations of the Richelot-Humbert construction do not apply in 
a straightforward fashion to the case of a genus 2 curve with complex roots. In the next 
section we describe the generalisation needed to apply this for our monopole curve. 

Consider the genus 2 curve C 

+ P{x)Q{x)R{x) = 0, 
P{x)^{x^a){x~a'), Q{x) = {x-b){x-b'), R{x) = {x - c){x - c'), 
where the real roots a, a', 6, 6', c, c' are ordered as 

a < a' < b < b' < c < c' . 
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We may associate to this triple of (real) polynomials {P,Q,R) another triple, {U,V,W), 
defined by 

(4.4) Uix)^[Q{x),Rix)], Vix)^[R{x),Pix)l W{x) ^ [P{x),Q{x)]. 

r 1 d/(x) , ^ dg{x) ^ ^ , ^ 

where [f,g\ := g[x) — f{x). The roots of the (quadratic) polynomials U,V,W 

are all real. If we set u < u' , v < v' , w < w' the roots of U, V, W respectively, then one finds 

(4.5) a<v<w<a'<b<w'<u<b'<c<u'<v'<c'. 

(Explicit expressions for the roots of U, V, W will be given below from which these inequalities 
can be proven.) Humbert gave a geometric perspective on this construction. Let Pi,Pi be 
the roots of (the quadratic) polynomial Pi (j = 1, 2, 3). We may view the six branch points 
{PijPi} G Eis six points on a conic Q. Now given a conic and six points lying on this we 
may construct six further points as follows. Consider the lines Li := Pip'i- The three lines 
Li form a triangle and the new points are the points of tangency to Q from the vertices of 
this triangle. These are the roots of [Pi,P2]. This is illustrated in Figure [s] below. 

We thus have a situation similar to Gauss' AGM case: to each pair of branchpoints one 
can associate another pair of points which are closer than the initial ones, and we expect a 
relation between the integrals of corresponding pairs. Iterating this process one shows for 
every pair the existence of a limit and obtains an expression for the integrals in terms of 
these limits. The relation between integrals suggested and proven by Humbert is 




Figure 8. Roots oiP,Q,R and U, V, W 
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and similarly for the integrals between the other pairs of branchpoints. Here A is the 
determinant of the matrix whose entries are the coefficients of P, Q, R in the basis (1, x, x'^), 
and S{x) is a polynomial of degree at most one. 

There is, however, an important element of difference between the elliptic and the hyper- 
elliptic cases. The map between elliptic curves (4.2) whose iteration leads to Theorem |4] is 



replaced by a correspondence in the hyperelliptic setting. A correspondence T : C — >■ C of 
degree d between two curves C and C associates to every point p € C a divisor T{p) of degree 
d in C, varying holomorphically with p |GH78| . A correspondence can be presented by its 
"curve of correspondence", Z — {{p,q) : q £ T{p)} C C x C' . In the case of the Humbert 
construction, the two curves are 

C: + P{x)Q{x)R{x) =0, C : Ay'^ + U{x')V{x')W{x') = 0. 

The correspondence between C and C considered by Humbert jHumOl] is of degree 2, and 
is given by the curve Z C C x C' of equations 



P{x)U{x') + Q{x)V{x')=0, 
yy' = P{x)U{x'){x-x'). 



(4.7) 

In analogy with the pull-back of a map one can also introduce for correspondences a 



linear map Sz ■ ^^{C) -^• ft^{C). Then eq. (4.6) can be interpreted as the relation between 
differentials 

(4.8) fe = Sild., 



together with an analysis of the image of the path joining a, a' (resp. 6,5' or c, c') under 
the correspondence. In fact this analysis, while mentioned in passing in |BM88j . is rather 
a crucial point in the extension of Richelot result to the case of complex conjugate roots. 
Even if one simply integrates on the right hand side along a straight line connecting two 



branchpoints, the image of this contour via the correspondence (4.7) may be considerably 



more complicated. Indeed one can obtain nontrivial homology cycles in the image and eq. 



(4.6) should be interpreted as an identity in Jac(C'), the Jacobian of C. One finds that 
certain half-periods in C are sent to periods in C' via the correspondence. 

With this background we may now state a version of the Arithmetic-Geometric Mean for 
genus 2 curves (with real branch points). 



4.2.1. The AGM method for genus 2 curves. Consider the genus 2 curve (4.3). Define six 
sequences (a„), (a^), (5„), (6^), (c„), (cj^) recursively by the conditions: 

• ao = a, Oq = a', bo = b, fej, = 5', cq = c, Cq = c'; 

• a„+i, a^+i, 6„+i, c„+i, are roots of C/„F„W„, ordered as follows 
(4.9) a„+i < < bn+i < 6^+1 < Cn+i < c'^+i, 

where, for every n, 

Pn{x) ^ {x - an){x - a'„), Qn{x) = (x - bn)ix - b'J, Rnix) = {x - Cn){x - c'J, 

Un{x)^[Qn{x),Rn{x)], K (a;) = [i?„ (a;) , P„ (x)] , Wn{x) ^ [Pn{x),Qn{x)]. 
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Bost and Mestre |BM88j give an explicit expression for these sequences: 



(4.10) 



with 



J'n+l 



"n+1 



Cn+1 



CnCn 


a^n^n 


~ Bn 


Cn + 


c'n ^ 0-n 


-a'n 


Kh'n 


- ana'n 


^Cn 


bn + 


b'n - an 


- a'n 


CnC-n 


- bnb'n - 


^An 



bnb'n 



^n+1 



51 





b'n - an - a; 


c c' 

'-n 


- bnK, - An 


c„ + 


c'n - bn - b'^ ' 


CnC-n 


— ana'n + Bn 



Cn 


+ C'n- 


bn - b'„ ' 


Wl+l — 






An = 


Vibn 


- c„)(6„ 


-<){b'n- 




- c' ) 


Bn = 


ViCn 


- an)(c„ 


- a'n) (c'n - 


-an){c'n 


-a'n), 


Cn = 


\/ian 


- 


-b'n){a'n 


- bn){a'n 


b'n)- 



These can be derived finding the roots for Un , Vn , Wn as foUows 



bnb'n T An 



(4.11) 



r -L r' — h — /)' ' 

^n ^ ^n ^n ^n 

CnC'n " QnO^ T -Bn 
Cn ^ '-n an ^n 
&n^n ^ flnfln "F C'n 



bn + - an - a'„ 

and ordering them according to (14. 91). One sees directly from the expressions above that 



Vn < Wn < w^ < Un < — ^n (^f- also cq. ( |4.5[ )). Thus we set 
(4.12) 



-^n+l 



Jn+l 



n+l 



■n+1 



^n+1 



and (4.101 then foUow. We then obtain: 

Theorem 5 (Richelot |Ric36j . Bost and Mestre |BM88p . With the above definitions, 
the sequences (an), (aJJ, (bn), (b'n), (c„), (c'^) converge pairwise to common limits 

lim ttn = hm a'^ — a = M(a, a'), 

n— Voo n— >-oo 

lim fe„ = Urn = = M(&, 6'), 
lim Cn — lim = 7 = Af(c, c'). 
Furthermore, for any polynomial S{x) of degree at most one, the following relations hold: 



(4.13) 



I{a, a) 
I{b,b') 
I{c,c') 



S{x)dx 



a ^-P{x)Q{x)R(x) 

S{x)dx 
b ^~P{x)Q{x)R{x) 
S{x)dx 
^-P{x)Q{x)R{x) 



irT 



= ttT 



= ttT 



S{a) 



[a - I3){a - -^)'' 

S{P) 
(/?-«)(/? -7)^ 

^(7) 



where 

00 

(4.14) T =[]<„, tn = 



{i-ol){i~ 13)'' 



n=0 



V (bn + b'n - On - a',^){Cn + c'^ - bn - b'J{Cn + < - a„ - fljj 
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The proof of the convergence of the sequences (a„), (aJJ (and hkewise (6„), (fe^); (c„), 
(c^)) is similar to that of the eUiptic case. Using (4.6) we find 



(4.15) 



S{x) 



dx 



S{x) 



t„ 



rdx 



■\/^ -Pn+ 1 Qn+1 -Rri 



=da;, 



and the relations (4.13) follow upon taking the limit for n — oo and using the residue 



theorem. Integrals between other pairs of branchpoints { e..g. a! and 6) may also be calculated 
using the same method in conjunction with appropriate fractional linear transformations. 
We remark that the integral between h and h' given above has opposite sign to that in 
|BM88| . because of a different choice of conventionf]^ 

4.3. Generalisation to the genus 2 case v^rith complex conjugate roots. We now 

generalise Richelot's method of the previous section to the case where the branchpoints are 
not all real but the polynomials P, Q, i? are still real. This corresponds to the three pairs of 
complex conjugate branchpoints, 



6' = 6, 



c, 



the case relevant for our monopole curve (1.5). We further order the roots such that 



Rc(a) = Re(a') < Re(6) ^ Re(6') < Re(c) = Rc(c') 
and for definiteness take Im(o) < 0, Im(6) < 0, Im(c) < 0. This splitting into complex 



conjugate pairs was given for the quotient monopole curve in (2.6) and (2.7). 



Of course all of the polynomial relations given in the previous sections extend to the case of 



arbitrary complex branchpoints and so the relation ( 4.8 ) between the differentials on C and C 



still holds true for this case. The difference with complex branch points arises at two points. 
First, with complex roots, there is no natural way to order the branchpoints and hence no 
natural way of splitting the branchpoints into pairs; thus there is choice in constructing a 
sequence of branchpoints to iterate. (This same feature is present with the ordinary AGM 
when the elliptic curve does not have real structure.) Second, as noted earlier, the image 
of the path between branchpoints under the correspondence may be quite complicated. 
Restricting attention to the case of the three quadratics P, Q, R having complex conjugate 
roots simplifies the problem somewhat. Although the initial branchpoints are complex and 



a relation analogous to (4.5) cannot be written, nevertheless the roots of C/q, Vb, VFq are real 



(as can be seen by considering their explicit expressions in (4.11)) and so can be ordered 



In contrast with the purely real case however, this ordering is not unique: in the real case 
the ordering of m, it', v, w', w, w' depended only on the relative ordering of a, a', b, b', 
c, c' on the real line; in the complex conjugate case this now depends on their imaginary 
parts. Depending on the ordering of it, u', v, v' , w, w' , equation (4.15) relating the integrals 



between the three pairs of branchpoints on C and C needs to be modified appropriately. We 
shall focus here on the case relevant for the monopole quotient curve. 



■^Bost and Mestre, in their note 2, p. 51 of |BM88) . claim that they want to recover the "classical 
identity" la — It Ic = 0. With our choice of convention for sheets, the relation between integrals become 
la + lb Ic = 0, which follows from the fact that the integral around a cycle encircling all the cuts, oriented 
so that the upper arc goes from negative to positive real values, is zero. 
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4.3.1. The AGM method for the quotient monopole curve. Consider the quotient monopole 
curve (1.5). Ordering the branchpoints as in section 2.1 we have 

(4.16) a = B4, a'^^a; b ^ B^, b' ^ B2; c = Be,c'^Bi. 



Calculating u, u' , v, w', w, w' via (4.111, and examining their relative ordering we find the 
following two cases: 

(4.17) case 1 : a > 0, v < w < w' < u < u < v'; 

(4.18) case 2 : a <0, u < v < w < u' < v' < w' . 



Case 1. This is exactly the same situation as considered in 



holds and so the sequences a„, aj^, 5„, 6^, c„, are still given by (4.10). Therefore the first 



BM88) . Thus (|4l2| stifi 



equality in (4.151 holds in view of (4.17) and so the integrals between complex conjugate 
pairs of branchpoints are still expressed by (4.13). We shall return to a discussion of integrals 
between other pairs of branchpoints shortly. 



Case 2. The different ordering of (4.18) means that (4.12) no longer holds for the first 
step of the recurrence. We modify this as follows. In view of (4.18) for n — 1 we take 

= V, bi — w b'l — 1 



fli = u, 



Cl 



w . 



Now at this stage the curve C has all real branchpoints and hence the Richelot-Humbert 
iteration can be applied as previously. Thus the only change in this case occurs at the first 
step of the recurrence. Let us denote by I{p, q) the integral on C between p, q on the first 
sheet, and by J^'^ the integrals on the curve C^*^ of equation + Pi{x)Qi{x)Ri{x) = on 
the same sheet. The integrals /^'^ can be then expressed by equations ( |4.13 ), using the 
AGM method for the curve C(i). We obtain an expression for the integrals I{p, q) on C as 
follows. Numericalljj^ we calculate the images on C under the correspondence (4.7) of the 
straight line contours of integration between the branchpoints of C The resulting contours 
of integration on the right hand side of (4.19) yield 

(4.19) 



I{a,a') = t^I {a\,bi), 

I{b, b') = io/'(c'i, ai) = t^{-l\a'^M) + Cl)), 
/(c,c') --to/'(fe'i,Ci). 



A similar numerical analysis of the images of the paths between other pairs of branchpoints 
then yields 

1 



(4.20) 



/(a,6)=-io (/ (ai,a'i)+/'(6'i,ci)), 
/(a',6')-^io (/'(ai,a'i)-/'(6'i,ci)), 
/(6, c) = (-/'(ai, a'l) + /'(a'l, &i) - I'{bM), 



lib', c') - -io (-/ (ai, a'l) " I'{a\M) - I\b,,b[)). 



"^More specifically, we study the images on the curve C' of equation Ay^ + U {x)V {x)W (x) = 0, in order 
to understand the first equality of | |4.15| l (as the second follows immediately from the first). We find, for 
instance, that the straight line between a and a' on C, call it 7(0, a'), is sent to a closed cycle encircling 
a'j and 61 on C' . Noticing that in the second equality of ( |4.15| l there is a factor of 1/2, absorbed in the 
definition of ti, we obtain that the image of 'y(a,a') on the curve + Pi{x)Qi(x)Ri{x) = 0, i.e. C' , is 
precisely the path from a'^ and bi. 
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Recalling that we are able to express the integrals I'{p, q) applying the AGM method to the 
curve C (with real branchpoints) using Theorem [5] and earlier remarks we are able then to 
calculate all integrals between branchpoints on C. 

Finally, we can apply similar numerical considerations to the integrals between non com- 
plex conjugate branchpoints in case 1 to obtain 



(4.21) 



/(a, a ) 


= hi (ai, a'l), 






I{b,b') 


- W'(6i,&'i), 






I{c,c') 


= tol (ci,Ci), 






I{a,b) 


= \to{I'{a\M)- 






I{a',b') 




-^'(ci,ci)), 




I{b,c) 


= i to (-/'(&!, 6'i 


-/'(6'i,ci) 


- I'(ci,ci)), 


I{b',d) 


^\to {l'{b,,b[) + 


/'(6'i,ci) + 


I'icuc'J)). 



5. Solving the Ercolani-Sinha constraints via the AGM 

We have shown how the Ercolani-Sinha constraints are reduced to finding the (a, g) 
{a/j3'^'^,-i/l3) such 

(5.1) 0=/^' Y^ = {X^ + aX + gf+A 



and for the cycle c given by (3.2 1 and (3.4 1 for a > and a < 0. Our strategy is as follows. 
Using the arc expansion (2.141 we may express the cycle c in terms of integrals between 
branchpoints. These integrals are then evaluated via the AGM method of the previous 
section using (4.13 4.211 for the case a > and (4.19 4.191 for a < 0. This has been 



implemented in Maple. The advantage of the AGM method is that it is much faster than 
the direct numerical integrations between branchpoints as it deals only with polynomial 
manipulations; moreover, the convergence of the sequences (a„), (aj^), (6n), (b'n), (c„) only 
needs very few steps, usually 6 or 7, for the precision we require. These considerations 
allow us to successfully solve the Ercolani-Sinha constraints iteratively, as described in the 
next subsection. A useful check of the method is the calculation of the matrix of a and 
b periods. These periods may again be reduced to integrals between branchpoints using 



the arc expansion (2.141 and consequently evaluated via the AGM. We find agreement with 
the same quantities being evaluated by other methods, yet with significant improvement in 
speed. 



5.1. Numerical solutions. To begin, we wish to find solutions to (5.1) starting from the 
point a — Q and g = 5^2. The integral over c is reduced to integrals between branchpoints 
as described above. We then proceed iteratively as follows. We start varying a by a small 
e, namely = i • e; we then vary 7 in smaller steps, 7^^^ = 5v^ + k ■ . For every such 
pair (ai, 7^^^) we calculate the periods using the AGM method, and hence compute the first 
constraint in (5.1). For every ai we take the ^i^k for which this constraint vanishes. We 
repeat this a sufficiently large number of times, obtaining the curve in Figure [21 We have 
used a step size of 10^^ for a S (0,2.8), and of 10^^ for a S (2.8,3.0) to obtain greater 
detail in this interval. For values of a e (0,3), the outcome is that we have a curve of 
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Figure 9. Solutions to the Ercolani-Sinha constraints 



solutions in the space of parameters passing through the points (0, 5y/2) up to the point 
(3,0). This point does not however belong to the solution curve. We may also extend this 
curve for negative values of a to the left of the point (0, 5-\/2). In Figurejojwe plot 100 points 
corresponding to a < 0, g > 0. 

The point (a, g) = (3, 0) is in fact a singular point as 4 of the branchpoints collide pairwise, 
giving two singular points at ±i. This results in a rational curve with equation 

(5.2) = {x^ + 3xf + 4 = (a;2 + A){x^ + if. 

Note that the curves of equation ( |5.1| with g ~ and a > are all hyperelliptic with the 
only exception being precisely the case a = 3 above where the curve is reducible. In terms 
of the original spectral curve (1.3 1 we have = 27/3^ and |HMM95] noted the loci under 
consideration here being asymptotic to this at one end. As a rational curve has no nontrivial 
cycles, the second of the Ercolani-Sinha constraints ( |3.6[ ) that fixes /3 tells us /3 = oo and so 
our solution curves are asymptotic to the point (a, g) — (3, 0). We remark that the behaviour 
of the solution curve in Figure [9] is consistent with the findings of [Sut97| , where Sutcliffe 
predicts that the curve ( |5.2[ ), describing a configurations of three unit-charge monopoles 
with dihedral D3 symmetry, constitutes an asymptotic state for a 3-monopole configuration 
(cf. eq. (4.16) in [Sut92]). 

When trying to extend the solution curve to g < 0, we do not observe any values of 
the parameters satisfying (5.11 with the first set of integers of (3.4). But since the point 
(a, g) = (3, 0) does not belong to the solution curve, continuity arguments do not prevent us 



using the second set of integers of (3.4). With these we manage to extend the solution curve 
through the point (0, —5\/2), which again corresponds to a tetrahedral monopole, now with 
a different orientation. This curve is also shown in Figure [9] We point out that the arc of the 
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-10 



Figure 10. p = /3(a) 



curve for g < is precisely the reflection with respect to the a-axis of the arc obtained for 
g > 0. Because of this symmetry we may focus attention on the case g > in what follows. 
Having determined the relationship between a and g we then use this to find /? via 



6/3 



1/3 



XdX 
Y ■ 



Again these are just integrals determined via the AGM and we present (3 — /3(a) in Figure 
[To) We shall interpret these results and compare them with other work in the next section. 



6. Discussion 

Before discussing our results let us summarise our argument thus far. In this paper we 
have constructed the spectral curve associated to charge three monopoles with cyclic (but 
no t dih edral) symmetry C3. By imposing cyclic symmetry the origi nal genus 4 spectral curve 
C ( 1.3 ) was shown to cover the genus two hyperelliptic curve C ( 1.5 ). By making use of a well 



adapted homology basis, the theta functions and data appropriate for the monopole solution 
were shown to be expressible in terms analogous data for the quotient curve C. Thus the 
construction of an appropriate spectral curve reduced to questions purely in terms of the 
curve C. In particular the transcendental constraints of the Hitchin construction reduce to 
the single transcendental constraint 
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for the scaled curve = {X'^ + aX + g)^ + 4 and a specified cycle c. This may be viewed 
as defining a function g = g{a) and the monopole curve is determined in terms of this. 
The "special function" g{a) warrants further study. Here we have made a numerical study 
of this. Our numerical study used a genus 2 extension of the arithmetic-geometric mean 
found by Richelot. Like its genus one counterpart Richelot's extension converges extremely 
rapidly and is an excellent means of evaluating such integrals. Richelot's method has (to our 
understanding) been used almost entirely in the setting of genus two curves with real branch 
points. When extended to the case of complex (here conjugate) branch points several new 
features arose. 

Our results extend work of both Hitchin, Manton and Murray |HMM95j and Sutcliffe 
|Sut97| which both describe cyclically symmetric charge three monopoles. In the former 
the following picture of the scattering of three monopoles, corresponding to geodesic motion 
along one of our loci is given. Three unit charge monopoles come in at the vertices of an 
equilateral triangle, moving towards its centre, in the xi — xi plane. Asymptotically this 
equilateral triangular configuration corresponds to the reducible spectral curve at a = 3. 
At a = the three monopoles coalesce instantaneously into a tetrahedron. Depending 
on whether the equilateral triangle is below or above the x\ — X2 plane we have distinct 
orientations of the tetrahedron corresponding to 5 > or g < 0. Finally the tetrahedron 
(with say g > 0) breaks up into a unit charge monopole moving along the positive x^-sxss, 
at (0, 0, 5) and an axisymmetric charge 2 monopole, moving along the negative xs-axis at 
(0, 0, —h/2). The reducible curve corresponding to the product of these configurations is 

= (r? + 2hO{v^ 2bvC + [b^ + y]C') = + - Sb'H' + 26(&' + ^X' 
from which we have the asymptotic behaviour at this end of the scattering given by 
(6.1) a- (7rV4-362), 7 - 26(62 + 7rV4), ^ _ 

where we have ignored terms vanishing as 6 tends to infinity. Sutcliffe investigated the same 
locus of monopoles numerically finding approximate twistor data by considering Painleve 
analysis of the Nahm data at the pole. This led to approximate forms of a, f3 and 7 described 
parametrically. In Figures [TTj [T2] wc plot this approximate data alongside the exact results. 
Despite not actually giving a spectral curve at any point the energy densities obtained by 
Sutcliffe upon solving the Nahm equations qualitatively refiect the scattering behaviour 



described above. To compare with the asymptotic prediction (6.1 ) we must revert to a and 
7. In Figure [13] we give a log- log plot of the exact values against this asymptotic prediction 
and alongside that of Sutcliffe's approximate curve. Our results reproduce the asymptotic 



behaviour of (6.II. In fact our approach could be extended to enable the calculation of both 
analytic and numerical corrections to this leading behaviour. 

There exists a rather nontrivial check of our results. We have argued that the Hitchin 
constraint H3 is automatically satisfied by o ur construction. This means that each of the 
"0 0] 



three theta functions 



- 
3 " 



{z;t) of (2.33) (with fc = 0, 1, 2) should be nonvanishing for 



A e (0, 2) and two of the three should vanish at the endpoints. A sample check is shown in 
Figure [T4| for a = —12.3, with Figure [TS] showing an enlarged portion of the k = curves to 
confirm the nonvanishing. Such a numerical plot is the only means we know of for verifying 
this condition. 

Although our construction is generically in terms of a genus two curve it may happen 
that this curve covers an elliptic curve. Such occurs for the tetrahedral monopole. Shaska 
and Volklein |SV04j describe when genus two curves give 2 : 1 covers of elliptic curves. For 
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Figure 11. g{a) for small 
values compared with Sut- 
cliffe's approximation (solid 
line). 



Figure 12. /3(o) com- 
pared with Sutcliffe's 
approximation (solid line). 




Figure 13. A log-log plot of the asymptotic behaviour of a versus 7 ac- 
cording to Hitchin, Manton and Murray (solid), Sutcliffe (dash) and here 
(dots). 

our curves these correspond to the axes a — Q (D3 symmetric monopoles) and 5 = and 
also to the solid line given in Figure |16[ We see that our curve covers an elliptic curve at 
two further points for < a < 3 and g > (and similarly for 5 < 0). These points do not 
appear to be otherwise special. For example, Sutcliffe studies the zeros of the Higgs field 
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Figure 14. A plot of the real and imaginary parts of the genus two theta 



functions (2.33) for fc = 0, 1, 2 




Figure 15. An enlargement of the fc = behaviour showing nonvanishing. 



and shows there is a point of positive a for which there is a 'zero anti-zero' creation event 
|Sut97j ; both our points are different from Sutchffe's. 
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